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The theoretical description of strongly correlated quantum systems out of equilibrium presents several chal- 
lenges and a number of open questions persist. In this paper we focus on nonlinear electronic transport through a 
quantum dot maintained at finite bias using a concept introduced by Hershfield [Phys. Rev. Lett. 70, 2134 (1993)] 
whereby one can express such nonequilibrium quantum impurity models in terms of the system's Lippmann- 
Schwinger operators. These scattering operators allow one to reformulate the nonequilibrium problem as an 
effective equilibrium problem associated with a modified Hamiltonian, thus facilitating the implementation of 
equilibrium many-body techniques. We provide an alternative derivation of the effective Hamiltonian of Hersh- 
field using the concept of an "open system". Furthermore, we demonstrate the equivalence between observables 
computed using the Schwinger-Keldysh framework and the effective equilibrium approach. For the study of 
transport, the non-equilibrium spectral function of the dot is identified as the quantity of principal interest and 
we derive general expressions for the current (the Meir-Wingreen formula) and the charge occupation of the dot. 
We introduce a finite temperature formalism which is used as a tool for computing real time Green's functions. 
In a companion paper we elucidate a generic scheme for perturbative calculations of interacting models, with 
particular reference to the Anderson model. 

PACS numbers: 72.10.Bg, 73.63.Kv, 72.10.Fk 



I. INTRODUCTION 

Quantum systems exhibiting an interplay of interactions 
and out-of-equilibrium effects are of significant interest and 
constitute an active area of research^— In contrast to the 
situation of equilibrium physics, however, there is currently 
no unifying theoretical framework for describing the dynam- 
ics of a generic quantum system out of equilibrium. The 
real-time Schwinger-Keldysh formalism^— has been suc- 
cessful in the treatment of specific systems^ 4 - but occasion- 
ally gives rise to pathological perturbative expansions which 
suffer from infrared divergences^ This can be viewed as 
an artifact of the infinite contour used in the integration, 
and one generally requires additional relaxation mechanisms 
to regulate the theory. While quantum impurity models 
in equilibrium have been extensively studied, 1 ^ probing and 
understanding properties of nonequilibrium steady states is 
a far more subtle task, many aspects of which are yet to 
be explored. Nevertheless, there have been significant ad- 
vances in our understanding via several distinct approaches, 
which include the scattering Bethe Ansatz^ 7 -^ field the- 
ory techniques^— time-dependent density matrix renormal- 
ization group (RG)r2iS time-dependent numerical RG j 23,24 
perturbative RG^r— Hamiltonian flow equations' func- 
tional RG'— strong-coupling expansions'^ diagram- 
matic Monte Carlo,— ~— and imaginary-time nonequilibrium 
quantum Monte Carlo42iiL 

Mesoscopic quantum objects such as quantum dots ("arti- 
ficial atoms"), characterized by a set of discrete energy lev- 
els, in contact with reservoir leads.r 2 . are described by quan- 
tum impurity models. In this context, the Anderson model, 
which describes a single discrete level coupled via tunneling 
to a Fermi-liquid sea, is of particular interest^ For a level 
with sufficiently low energy and with strong on-site interac- 



tion, the Anderson model mimics the situation of a Coulomb- 
blockaded quantum dot' In this regime, charge degrees of 
freedom on the level are frozen out and the Anderson model 
becomes intimately related to the Kondo model. The latter 
describes a magnetic impurity (manifested by the spin of the 
highest occupied level for an odd number of electrons on the 
quantum dot) entangled with the spins of a Fermi sea. (In 
the mesoscopic setting, the Fermi sea is replaced by the reser- 
voir leads.) The Kondo entanglement then produces a promi- 
nent Abrikosov-Suhl resonance in the density of states of the 
quantum dot at the Fermi level, which results in perfect trans- 
parency when the quantum dot is symmetrically coupled to 
its leads (the source and the drain)' 4 ^ 7 - Nanoscale systems 
can be routinely driven out of equilibrium by applying a bias 
voltage between the two reservoir leads. Among the various 
issues which arise out of equilibrium, the precise fate of the 
Abrikosov-Suhl resonance 4 ? when applying a finite bias volt- 
age, remains a delicate issue. In general, it seems essential to 
elaborate theoretical and numerical methods which will allow 
access to the full current-voltage characteristics of interacting 
mesoscopic systems. 

The present paper constitutes Part I of a two-part series. 
Its main purpose is to reformulate electronic nonequilibrium 
transport in quantum impurity models in terms of an effective 
equilibrium steady-state density matrix. This concept was ini- 
tiated by Hershfield'^ . who proposed that the appropriate 
density matrix in Boltzmann form could be constructed ex- 
plicitly by invoking the Lippmann-Schwinger operators of the 
system. Recently, there have been several numerical efforts at 
implementing this schem e 23 ' 24,40,41 . However, the theoretical 
foundations of this approach were yet to be firmly established. 
In this work, we concretely establish the validity of this effec- 
tive equilibrium approach by making use of the open system 
limit'which guarantees that the system relaxes to a steady 
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state. This well-controlled mathematical contrivance mimics 
the role of relaxation mechanisms, and does not require the 
explicit inclusion of the bath degrees of freedom at the level 
of the Hamiltonian. Furthermore, we demonstrate the equiv- 
alence between observables computed using the Schwinger- 
Keldysh framework and the effective equilibrium approach. 
It is important to emphasize the fact that the effective equi- 
librium description encompasses the case when one includes 
interactions on the dot. We then propose an imaginary-time 
formulation within this framework, which allows us to build 
general formulae and establish the mathematical machinery 
for evaluating quantities relevant to transport, and to sys- 
tematically implement perturbative and non-perturbative tech- 
niques familiar from finite temperature equilibrium field the- 
ory. In Part II,— we employ the formalism in the context of 
the nonequilibrium Anderson model and perform a systematic 
perturbative expansion in the repulsive interaction on the dot. 
Other aspects of the method, including numerica l appro aches, 
have been recently investigated in Refs. \2 lll23ll24ll40ll4 lll5 3ll54 
In addition, we note that the non-interacting resonant level 
model may also be a pertinent starting point in tackling sit- 
uations with strong interactions, as for example the large- N 
Anderson mode l 35 i 53 i 55 or the Toulouse limit of the Kondo 
models 

Our paper is organized as follows. In Section II, we in- 
troduce the model and define the "open-system limit" 
This concept provides a transparent framework that guaran- 
tees the existence and uniqueness of a steady-state density 
matrix and we will use it recurrently in our discussion. In 
Section III, we present the effective equilibrium approach and 
justify the validity of the method. We provide an alterna- 
tive derivation of the effective equilibrium density matrix pro- 
posed by Hershfield^ and establish its equivalence with the 
formal expression for the density matrix proposed by Doyon 
and Andrei 5 -^. Furthermore, we present useful recursion rela- 
tions for the expectation values of observables which underpin 
the equivalence between this approach and the time dependent 
prescription. In Section IV, we demonstrate the equivalence 
between observables computed within this framework and the 
Schwinger-Keldysh formalism which enables us to arrive at 
the Meir-Wingreen formula for the curren t 57 ' 58 . We then de- 
rive a compact expression for the charge occupation on the 
dot, whereby we identify the spectral function of the dot as the 
central quantity of interest. It is instructive to note that this ef- 
fective equilibrium formalism avoids the (often complex) cou- 
pled Dyson's equations for the various Green's functions on 
the Keldysh contour, and it is in principle possible to com- 
pute the Green's function of interest directly. In Section V we 
develop the methodology appropriate for the effective equi- 
librium approach in terms of a finite temperature imaginary- 
time formalism, which can be used to analytically compute 
real time Green's functions. 



II. MODEL 

As a prototype for our discussion of Hershfield's approach 
we consider a system consisting of two Fermi-liquid leads, 




FIG. 1: Schematic setup (upper panel) and energy diagram (lower 
panel) of a generic quantum impurity model out of equilibrium. The 
system is at a temperature T, and <& = /Lti — = eV denotes the 
voltage bias between the source and drain leads. The (bare) energy 
of the dot level is given by The energy broadening of this level, 
given by the width T, is due to the tunnel coupling between dot and 
leads. 



coupled by tunnel junctions to a central system with a number 
of discrete levels, the "dot", see Fig. Q] The Hamiltonian of 
such a generic system can be written in the form H = Hl + 
Hd + Ht- The first term 

H L = * a kJ aka c aka (1) 

aka 

describes the left and right leads (a = ±1), where c aka 
(fi ak ) annihilates (creates) an electron (strictly speaking a 
Fermi-liquid quasiparticle) in state k with spin projection a 
in lead a. The corresponding energy dispersion is denoted 
by e a k- These leads couple to the dot, whose Hamiltonian is 
given by 

a 

with d„ (d a ) annihilating (creating) an electron with spin a 
in the discrete level with energy e^. Any electron interactions 
are lumped into the contribution H mt , and we will assume that 
these interactions are localized on the dot (as it is the case for 
the Anderson model). Finally, the tunneling of electrons be- 
tween leads and dot is captured by the tunneling Hamiltonian 

Fr = -^E^( c i^+ h - c -)' o 

akcr 

where il is the lead volume (assumed identical for both leads) 
and t ak specifies the tunneling matrix element for electron 
transfer between state k in lead a and the discrete dot state. In 
the presence of a bias voltage, realized as a chemical potential 
difference $ = /ii — between left and right lead, the 
tunneling induces an electric current. 
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Typically, the steady state is reached after a short time de- 
termined by relevant relaxation rates in the leads. For the pur- 
pose of calculations aiming at steady-state quantities, it is con- 
venient to avoid the consideration of microscopic relaxation 
mechanisms and instead invoke the so-called "open-system 
limit" In short, this approach proceeds as follows: Ini- 
tially, up to some time t = to < in the early past, the tunnel- 
ing term is absent and the Hamiltonian of the system is given 
by Hq = Hl + Hu, such that leads and dot are decoupled. 
For t < t the system is hence described by the separable 
density matrix 



Po = exp 



-P [H 



(4) 



Here, j3 = (fcsT) 1 denotes the inverse temperature and 



J2ka ^aka Caka ^ S tne numrjer operator of electrons in 



lead a. Between times to < t < 0, the tunneling is then 
'switched on' adiabatically, i.e., H = Ho + HTe r,t 9(t — to), 
where the parameter 77 — > + defines a slow switch-on rate. 
At time t = the tunneling has reached its full strength, the 
system is in its steady state and observables can be evaluated. 

As demonstrated in Refs. HH and [TH me existence and 
uniqueness of steady state is tied to the validity of the inequal- 
ities vf/L <C \to\~ 1 "C T], where vf denotes the Fermi veloc- 
ity and L the linear system size. Intuitively, these inequalities 
ensure that hot electrons hopping onto a given lead at time 
to will not be reflected back and return to the junction be- 
fore the measurement process, and further that the process of 
switching on Ht remains adiabatic. We also note that the en- 
ergy scale of switch-on |<o| _1 suffices to smear out the energy 
level spacing vp/L. In this sense, the openness of the sys- 
tem provides the "dissipation" mechanism necessary for the 
steady state, allowing the high-energy electrons to escape to 
infinity and thus, effectively relax. On a more technical level, 
the inequalities result in the factorization of long-time corre- 
lation functions, which facilitates the proof of existence and 
uniqueness of steady state. 

The crucial ingredient in our discussion is that L/vf deter- 
mines the largest time scale in our problem. The exact proto- 
col by which we switch on the tunneling is irrelevant for the 
formation of steady state. Let us take as an example the non- 
interacting case, and for the sake of convenience assume that 
we have identical leads, i.e., t a k = t. Here, instead of adi- 
abatically turning on the tunneling suppose we do a quench, 
one observes that the transients decay with a relaxation time 
~ r, where T = 2nt 2 v denotes the linewidth of the dot. 59 
Here v symbolizes the density of states, which for simplicity 
is assumed to be a constant. To make the argument concrete, 
assume that the tunneling Hamiltonian is absent for t < 0, 
and the distribution functions in the leads are given by Fermi 
functions f(e — a%) (a = +1,-1 specify the source and 
drain reservoirs respectively). Furthermore, suppose that the 
dot is initially unoccupied, i.e. n™ = denotes the initial 
charge on the dot. The occupation of the dot at a time t(> 0) 



is given by 59 
n d (t) = n dss (l 



-2Tt\ 



2Te 



17 +/(*+«) 



cos [(u - e d )t] 
r 2 + (w - e d f 



Here 



n d ss 



du 



/(e-f) + /(e+f) 



r 2 



(5) 



(6) 



denotes the steady state expectation value of the dot occupa- 
tion, after the transients have decayed. It is interesting to note 
that if we turn the coupling of the QD to the left and right 
leads adiabatically to zero at the same rate, (formally this im- 
plies taking V — >• in Eq. ©) then the final QD occupation is 
given by 



fin 



f(e d 



<I> $ 

2 ) + /(£d + 2 ; 



(7) 



The fact that n™ ^ njj n , clearly illustrates the irreversibility 
of the turning on process. It encapsulates precisely how the 
system is driven out of equilibrium and is the reason we are 
required to analytically continue time to the Keldysh contour. 
This demonstrates that the steady state formalism cannot be 
continued from zero to finite tunneling perturbatively. There- 
fore, in our effective equilibrium approach, we include the 
contribution of the tunneling terms non-perturbatively. Inter- 
actions on the other hand can be turned on adiabatically. This 
procedure is reversible and does not exhibit similar anomalies. 



III. DERIVATION OF THE EFFECTIVE EQUILIBRIUM 
FORM OF THE STEADY-STATE DENSITY MATRIX 

The description of a nonequilibrium problem in terms of 
an effective equilibrium density matrix was first proposed by 
Hershfield in Ref. To facilitate such a description, the 
nonequilibrium steady-state density matrix p is rewritten in 
the usual Boltzmann form, 



p = cxp [-P{H - Y)} , 



(8) 



at the cost of introducing a correction operator, which we will 
call (following Hershfield's convention) the Y operator^ For- 
mally, the definition of this correction operator as 



Y 



H 



(9) 



is always possible. However, in this form it is neither 
particularly elucidating nor useful for calculating nonequi- 
librium transport properties. Hershfield put forward the 
idea that Y can be expressed explicitly and compactly in 
terms of Lippmann-Schwinger operator a 60 - 61 ipaka, which are 
fermionic operators that diagonalize the full Hamiltonian^ 



H 



(10) 
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Hershfield proposed that the Y operator has the general forrrpi is hence defined as 



(ID 



Since this operator encodes the entire $ dependence, the Y 
operator is also called bias operator. Note that Y vanishes at 
zero bias and the steady-state density matrix correctly simpli- 
fies to the equilibrium density matrix. 

In this Section we provide a detailed proof of the repre- 
sentation of Y given in Eq. (fTTT >. In contrast to Hershfield, 
we do not invoke the presence of additional relaxation mech- 
anisms to reach steady state. Instead, we make systematic 
use of the time-dependent open system approach.?^ ! which 
circumvents ill-defined expressions and makes the proof rig- 
orous. 



H - Y , n = 
= {h t -Y x , n = l 
-Y n , n>2 



(13) 



One can expand the exponential operator to collect terms in 
powers of the tunneling Ht such that 



1=0 ii,...,ij=0 



E 



^ Xii ■ ■ ■ Xii 



oo oo 

= EE 



J2 X H ...X tl =J2pn, d4) 



n— 1= Vii=n 



A. Proof of the explicit form of the Y operator 



where, following our general notation, p n denotes the order 
(Ht) 71 contribution to the steady-state density matrix. 



The structure of the proof is as follows. The starting point is 
the expansion of the steady-state density matrix into a power 
series in the tunneling Hamiltonian Ht- This can be accom- 
plished either by using the Hershfield form of the density ma- 
trix, or by employing the time-dependent framework of the 
open-system limit where the steady-state density matrix is ob- 
tained by adiabatically switching on the coupling in the far 
past. In both cases, one obtains analytical expressions for the 
power series in Ht- Order for order comparison of these ex- 
pansions then leads to a system of nested differential equa- 
tions for the bias operators. Finally, the solution to this system 
of differential equations is then shown to be identical with the 
representation of Y in terms of Lippmann-Schwinger opera- 
tors, see Eq. fTTY We should emphasize the fact that the ex- 
pansion in powers of Ht is a purely formal procedure, which 
we adopt for the sake of a systematic comparison, and the 
proof below is non-perturbative in Ht- 



1. Expansion of p using the effective equilibrium representation 

In the derivation of the explicit form of the Y operator, it is 
convenient to collect terms in orders of the tunneling Hamil- 
tonian Ht- We thus start by expanding Y = J2^Lo ^ n 
a series in powers of the tunneling, such that Y n oc (Ht) ■ 
In the following, the index n will always be used for power 
counting of the tunneling Hamiltonian Ht- From Eq. ([8]) we 
thus obtain 



p = exp 



exp 



n=2 



-/?|(i?o - Y ) + (H T - Yi) - J2 y ™ 

OO 



n=0 



(12) 



Here, we have regrouped the Hamiltonian with the Y operator 
order by order in the auxiliary operator X n ~ (Ht) 71 , which 



2. Expansion of p using the open-system approach 

Let us now employ the open-system approach (as outlined 
in Section II). In this case, the steady-state density matrix is 
obtained by switching on the tunneling in the early past to < 
0. Up to this time to, the leads are decoupled from the dot. 
The adiabatic switch-on of tunneling is facilitated by using 
HTe vt 0(t—to) for the tunneling Hamiltonian (note that in this 
scenario, the total Hamiltonian is therefore time dependent). 
At time t = 0, transients have decayed and time evolution 
has turned the original density matrix po = p(t = to) into 
the steady-state density matrix p = p(t = 0). We note that 
the condition of adiabaticity is strictly true only in the limit 
j^j <C Tj, where \to\ — > oo, which is assumed in the open- 
system limit.= We emphasize that the actual evaluation of this 
limit is deferred until the very end of all calculations. Keeping 
1/| to | and r) small but nonzero in the interim is crucial for 
mathematical clarity and for avoiding ill-defined expressions. 

This time, it will be convenient to work in the interaction 
picture (with respect to the tunneling), where in general 



Oi(t) = e iH o( f - t o)Q e - iH a(. t - t o) _ 



(15) 



denotes the interaction picture of the operator O. Note that 
the time to (not t = 0) has been chosen as the reference 
time where Schrodinger and interaction pictures agree, O = 
Oi(to)- [In this we differ from the conventions adopted 
by Hershfield^ who chose different reference times for the 
Heisenberg and interaction representation.] 

In the interaction picture, the density matrix satisfies the 
evolution equation 



d _ , . 
l dt~ Pl{t) 



[H T j(t),Pi(t)} 



(16) 
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This is formally solved by 

' (-»)♦* f n r rt 



Pi 



w = E^K n 



dti 



(17) 



x [Hr,j(ti),[fl"T,j(t2),| 

OO 



[Hr,/(*n),A)]] 



where the n = term simply fixes the boundary condition 
p(t = to) = Po and T denotes time-ordering of operators. 

It should be noted that Eq. (fTTT i already has the form of a 
power series in the tunneling and thus has been used to define 
the n-th order contribution p n ,i(t) to the interaction-picture 
density matrix. For later purposes it is useful to note that the 
contributions can alternatively be obtained from 



dt 



Pn,l{t) = -i[H TJ (t), p n -u(t)}, 



(18) 



i.e., a system of nested differential equations associated with 
the boundary conditions Po i(to) = Po an d Pn,i(to) = for 
n > 1. 



3. Derivation of differential equations for Y„ 

For the results to be consistent, we require that the steady- 
state density matrix p be identical to the steady-state density 

I 



matrix obtained in the open-system limit. This allows one to 
determine the correct form of the Y operator, order for order 
in the tunneling. 

We thus impose the identity of the steady-state density ma- 
trices 



JH {t-t ) pne -iH {t-t ) = p nI {t), 



(19) 



where we have transformed p n into the interaction picture. 
Eq. ( fT9l is expected to hold for all t > 0, since at that 
point the interaction has been fully switched on, and the time- 
independent Hamiltonian (used in Hershfield's effective equi- 
librium approach) and the time-dependent Hamiltonian (from 
the open-system approach) are identical. 

We now differentiate the left-hand side of Eq. (fT9b with re- 
spect to time and use Eq. ( [Pit to obtain 



dt PnAt) 



(20) 



E 

i=i 



E E^ 

iiH Hi=n fe=l 



dX lkJ (t) 



I J t V,.;. 



where X ih j = e iH °^~ to) X ik e~ iHo( - t ~ t °'> denotes the interac- 
tion picture representation of Xi k . Similarly, we may differen- 
tiate the right-hand side of Eq. ( fT9l . The resulting commuta- 
tor is given in Eq.([T8l and contains p n _ 1I {t)^ for which we 
substitute the corresponding expression from Eq. ( fT4b . This 
way, we obtain 



-Pnj{t) - i[ Pn - h M H T j(t)} = lJ2 



dt 



OO 

E 

(=1 



J2 E x ^i(t) ■ ■ ■ l x iM H T,i(t)} ■ ■ • x nJ (t) 

I 

^X il , I {t)...[X ih . ltI (t) t H TJ (t)]...X il , I (t), (21) 



ii+«2+— +i(=n k=l 
I 



where for n < we define X n = 0. In the last step of Eq. 
( f2Tb . the summation constraint is shifted from n — 1 to n to fa- 
cilitate the comparison with Eq. ( f2Qb - This comparison yields 
the relation 



valid for any operator O, the boundary condition p(t = to) = 
po implies^ 



lim Kj(t) = - Va7V a . 



(24) 



-X n!l = i[X n _ ltI (t),H T!l (t)]. 
dt 



(22) 



It now remains to prove that the following interaction-picture 
Finally, utilizing the relation ( TOT l between X and Y operators, expression of 



one finds that the Y n operators also satisfy 
d 



™ =|Z)«V'La,/(*)V'afca I /(t), (25) 



dt 



Y nJ = i[Y n _ ltI (t),H T!l {t)}. 



(23) 



in terms of the Lippmann-Schwinger operators i^ aklJ j (t) 
To obtain the Y operator from these differential equations, it represents the solution to the above initial-value problem, 
is crucial to specify the boundary conditions at the initial time Once we have recapitulated the crucial properties of these 



t = t Q . Given the relation e iH {t-t ) e O e -iH {t-t ) 



Lippmann-Schwinger operators in the following subsection, 
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it will be simple to confirm that this ansatz indeed solves the 
differential equation ( |23l subject to the boundary condition 



4. Properties of the Lippmann-Schwinger operators 

It has been shown^ for the generic model defined in Sec- 
tion II that the Lippmann-Schwinger operators are fermionic, 



{V'ifecri V'a'fe'cr'} = 8 aa >6kk'Saa', 



(26) 



and diagonalize the full Hamiltonian (including all interac- 
tions), 



(27) 



One can expand ip ak „ in powers of the tunneling Hamiltonian 
Ht, i.e., 



4'\ = v ^ t 



(28) 



n=0 



such that Vdfccr „ oc (Ht)"- In the interaction representation, 
the operators tp^u- satisfy 

= ilH^^jit)} - i[H TJ (t),^ aklTiI (t)], (29) 

subject to the boundary condition 

lim ipak<T,i(t) = c aka . (30) 
t\t 

Equation d29l can be further simplified, and collecting orders 
of (Ht)™ 1 , cast into the set of nested differential equations 



dt 



1-1 (31) 

As discussed in Ref. [53l the formal solution can be ex- 
pressed compactly as 



r OL 



1 



ka aka 



e ak - C + ir] 



TC a kc 



(32) 



in terms of Liouvillian superoperators C and Ct- The ac- 
tion of such superoperators on any operator O is defined as 
CaO — [A,0]. Further, rj — >• + is a regularization similar 
in spirit to the one utilized in the more familiar Lippmann- 
Schwinger states in scattering theory. From Eq. ( l32l one can 
write the detailed form of the operator ^L,_ n 



^Ika = 



n=0 



1 



e a k - Co + ir\ 



-Cn 



~aka 



— ^aka, Ti- 
ll 

(33) 



Note that when the tunneling is set to zero, one finds indeed 



5. Conclusion of the proof 

It is now simple to verify that the ansatz for Y n j(t) solves 
the set of differential equations given by Eq. d25l) . and obeys 
the appropriate boundary condition d24t : We start by writing 
the (HtY 1 contribution to the Y operator as 

<|> n 

Y n j(t) =^a-J2 i>ik*, P j(t)^k*,n- P j(t). (34) 

aka p—0 

Differentiating this with respect to time yields 



d_ 

di' 



Y nJ (t) = J2 (r^ikavJ (*) J i>*ka,n-p,l(t) 

aka p=0 L V ' 



+ V'La,p,j(*) ( ^aka,n-p,l{t) 



, (35) 



so that in conjunction with Eq. d29l one obtains 

i j*. n—l 

d v— \ $ v — + 



dt 



aka p=0 

= i[Y n _ 1>I (t),H t ,i}. 



(36) 



This concludes the proof. 

In summary, we find that the steady state dynamics of the 
system can be described by an effective equilibrium density 
matrix of the form 



where the operator 



p = e 



Y =jJ2 ^aka^ka 



(H-Y) 



(37) 



(38) 



encodes the entire nonequilibrium boundary condition of the 
system. 



B. Recursion relations for expectation values of observables 

In addition to the previous proof, we follow Hershfield^ 
and underpin the equivalence of the adiabatic approach and 
the effective equilibrium approach by showing that they lead 
to identical recursion relations for expectation values of ob- 
servables. Again, the systematic use of the open-system limit 
makes the proof sound. 

Let O denote a generic observable. To derive the first recur- 
sion relation, we will decompose the steady-state expectation 
value (O) into a series counting the powers of the tunneling 
Hamiltonian Ht- In the first step, we thus substitute the ex- 
pansion ( TBI of p, 



(O) 



tv\pO] _ trlZZoPnO] 



(39) 
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Pulling out a factor of l/tr[po], expanding the denominator 
as a geometric series, and collecting terms order for order in 
Ht, we can rewrite this as 



<0> =££(-*)' £' n 

.,j ( =ls=l L 



n=0 ;=o 
oo 



.a- 



tr[po 



tr[Po] 



(40) 



where denotes a restricted summation, subject to the con- 
dition Y^ s =i Js < n. Eq. (1451 allows one to prove the impor- 
tant recursion relation 



(0)n = 



tr[po] 



(41) 



which relates the 77-th order term to the expectation value of 
with respect to p n , and all lower-order terms (0) m (m = 
Q,...,n- 1). 

Now, we turn to the time-dependent representation assum- 
ing an adiabatic switch-on of the tunneling. The steady-state 
expectation value of an observable may now be obtained via 



(0) 



tr[fr(0)O/(0)] 



(42) 



where we have switched to the interaction picture. To arrive at 
the desired recursion relation, we again decompose this into a 
power series of the tunneling Hamiltonian, 



tr[po] 



(43) 



n=0 



where 



(O)n 



H)« 1 



tr 



" r rt=0 



77! tr po 
x [H T ,i{ti),[H T .i{t2 



(It, 



(44) 



[ffr,/(*n),A)]]---]O/(0) 



This expression for (0) n allows one to derive the second re- 
cursion relation, which reads 



{(D)n = tr[p n j (0)0/(0)] _ y, tryo| (0>n _ ( 



k=l 



(45) 



The details of this derivation are given in Appendix [A] The 
proof makes explicit use of the factorization of two-time cor- 
relation functions in the limit of large time separation^ The 
agreement between Eqs. (l4Tb and (l45l l underpins the equiva- 
lence between the time-dependent adiabatic approach and the 
effective equilibrium approach. 



IV. CORRESPONDENCE WITH THE 
SCHWINGER-KELDYSH APPROACH 

Having obtained the effective equilibrium form of the 
steady state density matrix, we now illustrate how one can 



compute transport observables such as the current and the 
spectral function within this description. We propose an imag- 
inary time formulation for treating such nonequilibrium sys- 
tems in a manner similar to finite temperature field theory. 

In imaginary time, we define the propagation of an operator 
by 

0(t) = e^ H - Y ^Oe-< H - y ^ = e T ^ c ^O. (46) 

The nonequilibrium thermal Green's function is defined on 

< r < f3 as 

G 0l o 2 (r) = -<T[0i(t)0 2 (O)]) = -<Oi(t)O 2 (0)>. (47) 
Fourier transforming in imaginary time this results in 



iuj n — C + Cy 



02 



(48) 



where cj n = (2n + l)ir/(3 {11 £ Z) denotes the fermionic 
Matsubara frequencies. 

Switching to real time, the Heisenberg representation of an 



operator O is given by 0{t) 



AHt 



Oe 



-iHt 



The nonequilib- 



rium real time retarded Green's function can then be expressed 



G% l02 (t) = -i9(t)({Oi(t),O 2 (0)}) 



-i0{t) 



tr [ e -^ H ~ Y ) {0x^,02(0)}] 
tr [e-K B -*)] 



(49) 



By using the spectral representation and then Fourier trans- 
forming, we obtain from this 



C + 47/ 



(50) 



We emphasize that, despite the effective equilibrium character 
of the Hershfield approach, there does remain one important 
difference between the effective equilibrium description and 
any regular equilibrium many-body theory. This distinction 
arises from the different propagators in imaginary versus real 
time; namely, in Hershfield's effective equilibrium formalism 
all Heisenberg operators in imaginary time evolve under the 
modified Hamiltonian H — Y, whereas Heisenberg operators 
in real time evolve under the Hamiltonian H. As a result, 
imaginary time and real time Green's function, GoxO^ii^n) 
and Gq (id) are not simply related by a direct analytic con- 
tinuation i(d n — >• id + ir\. 

As an observable of prime interest in transport, let us now 
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consider the current flowing through the dot. It is given by 

h + i-i e/c^(iVl(^)-iv_l(^)) , 



i = 



2 

dt 



dt 



aka 



aka 



'2V?2 



= Im 



.afccr 



it(T = 0) 



(51) 



where — e denotes the charge of the electron. Here, we have 
used the fact that the system is in a steady state and the cur- 
rent is unchanged under time translations. As we shall see 
below, it will prove to be most convenient to express the cur- 
rent in terms of the Fourier representation of the imaginary 
time Green's function, namely 



Im 



Eetak 1 r , . , 

.a.kauj n 



(52) 



We now show the equivalence of this approach with the 
Schwinger-Keldysh formalism and recover the familiar Meir- 
Wingreen formula for the steady-state currenti^iS For sim- 



plicity, let us assume that the coupling to the dot is inde- 
pendent of k, i.e. t ak = t a , and that the leads are identical 
fat = efc. (This is not a strict requirement and the proof can 
be easily generalized.) Our starting point is Eq. d52l and we 
now write explicitly 




Using Eq. ( l32l for c aka and evaluating the effect of Ct, we 
obtain 



4>aka, 



iuj n — C + C\ 



-dt 



e k + C — ir] iu> n — C + Ly 



-dl 



(54) 



With a small trick, one can show that the second term does 
not need to be evaluated when computing the current I. To 
see this, we exploit the fact that in steady state, the current 
in the left and right junctions must be identical and hence, 
we can write the current as a weighted average of the form 



I = (t?I_i +t 2 _ 1 I 1 ) I (t? + fjjj. This eliminates the pres 



ence of the second term in Eq. (l54l from the expression for 
the current. Proceeding with the remaining term we obtain in 
several steps, 



J 




Ipaka, 



iu n — C + Cy 



dl 



1 



E 



iu> n + C — Cy 



1pak<r,dl 



e k + C- 

r 



where the transition from the first to the second line is facil- 
itated by the general relation (1B21 i. Also, in going from the 
sum over k to the integral over e k we make the replacement 
ftJ2k~^ I v d- e k- Here we have followed the standard proce- 
dure and assumed the spectrum has been linearized. Further 
v denotes the density of states, which is assumed to be a con- 
stant. Finally, this results in the well-known Meir-Wingreen 
expression for the currentjSS 



/ =2e- 



rir_! 
Ti + r_ x 



de k A d (e k ) 



- f Ufc - 



(56) 



where 



A d (e k ) = --Vim 



(57) 



is the nonequilibrium spectral function of the dot and is gen- 
eral a function of the bias voltage $. We shall however sup- 
press this explicit bias dependence of the spectral function. In 
the above expression T a = irt^y denotes the partial broaden- 
ing of the level due to the coupling to the lead a. 

Let us now derive the above identities via a slightly dif- 
ferent route and recover the traditional form of the Meir 
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Wingreen formula. Proceeding as in Eq. ( l54b one observes part of the above expression 



iuj n — C + Cy 



-dl 



txk 



(58) 



\ I iuj n + C — Cy e k — C + ir/ 
It is simple to show, in a sense similar to Eq. d55l l that the 1st 



Thus, using Eq. ( f5TT > together with Eqs. d55l l and d59l ) we get 



' = «X«^(^et to (r = 0)-^ 4 (r = 0) 



2tt 



/ <fc* / e fc - a- Gf 4 (e*) - G£ 4 (e fe ) 



-dl 



o iw„0 + 



iuj n + C — Cy " , e/c — C + ir/ a \ I \) €k + C — ir] a ' iu> n — C + Cy a 
Using an identity from AppendixlBl while integrating over e k in the last 2 terms, we find that 



/ = ^E« r « 



de k f \e k - a- 



4> 



(60) 



(61) 



Now 2inQ d a (t = 0) = 2in (d\d a } = 2nG < (t = 0) = de k G < (e k ). This helps us recover the traditional form of the 
Meir Wingreen formula: 



I = i— aT r 



de k < f e k - a 



(62) 



r 



We proceed in a similar fashion to derive an expression for 
the charge occupation on the dot 



E<4^> = ^E^4(* 



de k A d (e k ) 
de k A d (e k )f s (e k ,$), 



-§)+/(*+! 



where 



1 



/ ( <* - t ) + / ( £ fe + 



(63) 



(64) 



Note that on taking the limit $ — > 0, Eq.d63l reduces to a 
form well-known from equilibrium many -body theory. In Ap- 
pendix[B]we provide a straightforward derivation of this result 
using the spectral representation within the effective equilib- 
rium formulation. This result can also be obtained within the 
Keldysh framework in a slightly indirect manner as we have 
shown in Appendix [C] This perhaps accounts for its omis- 
sion in the existing literature, where the charge occupation is 
typically obtained using the lesser Green's function G < . 



We thus identify the non-equilibrium spectral function 
(w, $) as the central quantity of interest, which one can use 
to calculate transport observables. In the effective equilib- 
rium scheme this can be done directly, without resorting to 
the Dyson equations which couple the various Green's func- 
tions, and can be complicated. It is instructive to note that 
for the non-interacting case the spectral function is indepen- 
dent of the bias, but the introduction of interactions generally 
imparts to the spectral function a complex bias dependence. 



V. GENERATING FUNCTIONAL FRAMEWORK 

As usual, the evolution in imaginary time of an oper- 
ator in Heisenberg picture is governed by the propagator 
e~ T ^ H ~ Y \ which closely matches the form of the density ma- 
trix e~P( H ~ Y '. This similarity makes the imaginary time de- 
scription very convenient for the evaluation of Green's func- 
tions when using functional techniques. We will now outline 
the general methodology for constructing Green's functions 
using the effective equilibrium density matrix within a func- 
tional framework. 

We construct the generating functional as a coherent state 
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functional integral over Grassmann variables, corresponding 
to the Lippmann-Schwinger operators 



z[j,j*]= f i[(vr ak ^j a 



ktT) e- s V> J '\ (65) 



where S denotes the action 



S 



ah<7 



i'lkai^r + take, - a^/2)lp aka . 



(66) 



and for notational convenience, we have suppressed the fact 
that the Grassmann variables are functions of imaginary time. 
Any Green's function in terms of the Lippmann-Schwinger 
operators is now given by an appropriate functional derivative 
of the generating functional 



(r) 
1 



-<T[Vw(r)^ fev ,(0)]) 
5 2 

Z(J,J*] 



Z 5J* k(7 (T)5J a , k , a ,(0) 



(67) 



j=o 



Here Zo = Z(Q) denotes the generating functional with the 
source terms absent. 

We observe that the generating functional and consequently 
the Green's functions have a trivial form in terms of the 
Lippmann-Schwinger operators. However for transport we 
have to compute Green's functions in terms of the original 
degrees of freedom of the system such as Q , a and Q a . 

&ct@'ct Cakcr O'er 

The goal now is to express the c aka and the d a in terms of 
these Lippmann-Schwinger operators. This is by no means a 
trivial task, since the precise form of the Lippmann-Schwinger 
operators is in general not known. However for a system 
governed by single particle dynamics, the exact form of the 
scattering states is readily calculated and the actual degrees 
of freedom of the system and these states bear a linear rela- 
tionship to each other. In this limit, it is straightforward to 
pass from one set of operators to the other, and it will consti- 
tute the starting point of our calculations. Later this will lay 
the foundation for computing Green's functions in interacting 
theories, perturbatively in the interaction. 

For systems governed by single particle dynamics the 
Lippmann-Schwinger operators are given as a linear combi- 
nation of lead electron and dot electron operators 

^aka= E A «'^' C l' fe V'+E K -'4': (68) 
a'k'a' a' 

where A a 'k><r< and k„i are appropriate c-numbers. 

For interacting theories, the form of the Lippmann- 
Schwinger operator turns out to be quite non-trivial and can 
be written schematically as 

i'l = E A w...«i...€i...ft- 

Vi ,£j ,v' k ,f ( 



where the summation is over v i, £j, v' k , £| € {0, 1} with 
the constraint £\ v t + £\ £3 - X)k v' k - J2 t = 1, so that in 
total there is one more creation operator than there are anni- 
hilation operators. Here the coefficients A Vl ___ v i ...^ ...f... are 
c-numbers and we have used an abbreviated notation i,j,k,l 
for the labels of the single-particle lead and dot states. Thus, in 
a generic interacting model a simple way to pass between the 
Lippmann-Schwinger basis and the original basis is in general 
not easily accomplished. The construction of the Lippmann- 
Schwinger operator order by order in the interaction is how- 
ever well defined and from this we can compute the Green's 
functions in powers of the interaction strength, as shown in 
the companion paper (part II). For the purposes of numer- 
ics, it is possible to circumvent explicit construction of the 
Lippmann-Schwinger operator by working at the level of the 
density matrix^ 

The functional formalism described above provides a con- 
venient platform for the implementation of non-perturbative 
techniques. We are currently using this scheme for a large- 
TV treatment of the infinite-?/ Anderson model out of equi- 
librium, following the slave-boson language of Read and 
Newn o 35 ' 55 . The saddle point solution of the functional in- 
tegral in this case, mimics the non-interacting theory and it is 
possible to systematically build corrections in powers of 1 /N. 



VI. CONCLUSIONS 

In this paper, we have provided a rigorous alternative 
derivation of the effective equilibrium density matrix ap- 
proach proposed by Hershfield, in the context of nonequilib- 
rium quantum transport using the concept of an "open sys- 
tem". Furthermore, we have illustrated the methodology of 
this formulation for a quantum dot, in which the reservoir 
leads play the role of good (infinite) thermal baths, such that a 
unique steady state exists. The latter can be cast into the form 
of an effective equilibrium density matrix, where the associ- 
ated modified Hamiltonian can be explicitly written in terms 
of the Lippmann-Schwinger scattering state operators. 

We have detailed the foundations of the theory and demon- 
strated a rigorous correspondence between the observables 
computed using the Schwinger-Keldysh formalism and the ef- 
fective equilibrium approach. We emphasize that an advan- 
tage of the effective equilibrium approach is that it is con- 
ducive to the implementation of numerical methods, such 
as numerical R G 23 24 and Quantum Monte Carlo. 41 Further- 
more, we have built a generating functional framework using 
an imaginary-time formulation that will constitute the basic 
mathematical machinery for our computation of the steady- 
state dynamics of the system. In particular^ 2 - we will de- 
velop an infrared-convergent systematic perturbation theory 
scheme to treat interactions within the effective equilibrium 
approach. This scheme enables a direct computation of the 
non-equilibrium spectral function, from which we can in- 
fer transport observables, without resorting to (often compli- 
cated) coupled Dyson's equations which arise in the Keldysh 
context. 

Finally, we underline that the non-interacting resonant level 
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model may serve as an appropriate starting point in the 
case with strong interactions, e.g., in the large- N Anderson 
mode l 35 ' 55 or in the close vicinity of the Toulouse limit in 
the Kondo model. 5 * The effective equilibrium method may 
also be generalized to more sophisticated systems, such as 
one-dimensional mesoscopic systems embodied by the Lut- 
tinger paradigm^ 5 - where it is possible, in principle, to com- 
bine bosonization and the scattering pictured 
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Appendix A: Derivation of the recursive form of O r > 



We can cycle p p j(to) in the above expression out of the 
nested commutators, thereby inverting the nested commuta- 
tor structure to re-express this as 



tr p P ,i(t )[. . . [0,(0), H T ,i(t n )}, ...H TJ (t p+1 )} . (A5) 



In this appendix we provide the detailed derivation of 
the recursion relation d45l l when working within the time- 
dependent approach, switching on the tunneling adiabatically. 
Our starting point for the derivation is the expression for the 
order-(iJy) contribution to the steady-state expectation value 
of some operator O, 



(o) n =L*5ta 



dtr, 



dt n -i 



dti 



tr lPoJ Ut J to J t 

[H T j(t n ), [fr Ti/ (t n _i), [. . . [HTjih), po}} . . .}0 T (t = 0) 

(Al) 

see Eq. d44b in the main text. Carrying out the t\ integration 
and using Eq. ( fToT l. one obtains 

<0> " = W" 



/•*„ rt 3 

dt n I dt n -! ■ J dt 2 (A2) 

to J to J to 



x [H TJ (t n ), [. . . [flr,i(*a). (Pi,/(*a) - Pi,i(h))} ■ ■ .]Oi(0) 

We now retain the term involving pu(to), and proceed by 
carrying out the ti integration for the term containing pij (£2). 
This integration yields 2 terms: one with p2,i(to) which we 
retain, and the other with P2,iitz) which we subject to further 
integration. Continuing in this fashion until all the variables 
ti , . . . , t n in the latter term have been integrated out, we arrive 
at the result 

1x^/(0)0/(0)] 



(O)n = 



tr[/?o] 
(_j\n—p /•tn+i=0 

to 



(A3) 



dt, 



In 



dt, 



dt 



x tr I [flr, j(t„), [-.., [H T j(t p+ i),P P ,i(t )} . . .]Oj(0)| 



(A4) 



This allows us to rewrite (0) n in the form 



(O)r 



tr \p n ,i(0)Oi(0)] 
tr[po] 



(A6) 



V \ ' / dt n / dtn-1 ■ 
tr tP0] J t0 J t 



dt p +i 



x tr 



p P ,i(t )[. . . [Oi(0),H T ,i(tn)], ■ • ■ , H T ,i(t P+ i)} 



To proceed further, let 



«••■»(> = 



tr ( e -/3(H -Fo)...) 

tr ( e -/3(ffo-i-o)) ' 



(A7) 



denote an expectation value with respect to the density ma- 
trix po . It has been shown by Doyon and Andrei in Ref . EU 
that, in the open system limit, long-time correlation functions 
factorize, i.e. 



«O(ti)O(*2))) -> ((O(*i)»o«0(t 2 ))) (A8) 



as \t\ — t%\ —> 00. We will now employ this factorization 
in Eq. JA61 >. Renaming the integration variables t\ = 
t' 2 = t p+2 , . ■ ., t' n _ p ee t n and inserting p Q p^ x = 1 we get 
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(O)r 



tr [p„,/(0)O/(0)] 



(A9) 



(_j)«-P 



E 



dt' 



dti 



n—p—l 



dt[ tr \ Po^PpM [[• ■ ■ [O/(0),gr,j(4- P )], ■ ■ Pt,;(<'i)] 

^(*o) B(t' lt ...,? n _ p ,t=0) 



We now recall that in the open system limit we take vp / L <C l/|io| "C t? with to - ^ — oo. This implies that the e vt term in 
i?T,/ essentially cuts off the lower limit of the time integrals at a time ~ \Jr\ much earlier than to. Thus, Eq. dAlOt can be cast 
into an integral over ((A(t )B(t[, . . . ,t' n _ p ,t = 0))) , where \t[ - t Q \,. . .,\t' n _ p - t \,\0 - t \ -> oo. Using Eq. dA8l ) we get 



and hence 

(O)n 



({A(t )B(t[, . . . , t' n _ p , t = 0))) = {(A(t ))) {(B(t[, . . . , t' n _ p , t = 0))) , 
tr [ PnJ (0)0/(0)] 



tr[po] 



(A10) 



(All) 



j2 tejppjito)} tI 



P =i 



(-0 



•sn-p 



n—p 



dt. 



71— p— 1 



to 



<[[. . . [O;(0),%(C,)] ■ ■ • ,}HT,l(t'l)}P0 



Inverting the nested commutator structure as previously, we get 

tr [/w (0)0,(0)] 



(O)n = 



tr[p ] 
tr[p p ,j(f )] ^ 



P =i 



tr[po] 



dt! 



to 



n—p 



dti 



to 



n—p—l 



( i(' 1 [%((„l[...,[%(i' 1 ),p„]] 



(A12) 
10,(0) . 



to 



Finally, comparing the second term in the latter equation with 
Eq. (IA1 b and noting that tr[p p j(to)] = tr[/j p /(())], we arrive 
at the final recursion relation 



(O)r 



tr [pn,j(0)Oj(0)] 
tr[po] 



tr[ Pp , 7 (0)] 
^ tr[po] 



(O)n-p- 



(A13) 



Appendix B: Derivation of n d using the spectral representation 

In this appendix we provide additional details on the equiv- 
alence between the effective equilibrium approach and the 
Schwinger-Keldysh formalism and explicitly show the iden- 
tity of the electron occupancy on the quantum dot as obtained 
in the Schwinger-Keldysh formalism and as obtained with the 
effective equilibrium approach. We start by stating two useful 
identities involving Liouvillian superoperators. The first iden- 
tity regards the imaginary-time Heisenberg representation of 
an operator O and is given by 



0(t) = e T(H - Y ^Oe- T(H - Y ^ = e T ^ c - c ^O. 



(Bl) 



The second identity is a small trick which transfers a Liouvil- 
lian superoperator from one side of an anticommutator to the 



other, 



z-C 



:0, 



z + C 



o u o 2 



(B2) 



where z is a c-number. 

We now set out to demonstrate the equivalence of the 
Schwinger-Keldysh expression for the dot occupancy, 



nd 



1 f°° 

-J2 deA da (e) 

Moj^ AcUefc) 



/K + t) + 'K 

* + 1) + / U - 1 



(B3) 



and the expression obtained in the effective equilibrium ap- 
proach, 



£<4^> = £^ 4 (t = o). 



(B4) 



It is crucial to recall the definition of the Green's func- 
tions in the effective equilibrium approach (see Section HVb . 
The Fourier representation of the (imaginary) time ordered 
Green's function is given by 



-Oo 



(B5) 
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and the retarded/advanced real time Green's function in fre- where, for simplicity, we are considering symmetric tunnel 
quency space is obtained by couplings and have set t\ = t_i = t. Similarly, one finds for 



^tl(«) = ({Ou Z ^-0 2 }). (B6) 



the advanced Green's function 



To prove the equivalence between Eqs. ( IB3I ) and ( IB4b . we 
start with the Keldysh expression and evaluate the spectral 

(B9) 



function 

1 i r 

^radv 



Utilizing these relations in the evaluation of the occupancy, 
we find 



(B7) 

With Eq. JB6b we evaluate the retarded and advanced Green's rid = (BIO) 



functions involved, 



1 1 / , r ,+ , , , . X „ / $ 



^(({4>^})-h.c.)/ 



(B8) Writing the Fermi function as a Matsubara sum, one finds 



J 



E f ( 1 4. ,-,„„ X r _ r3 ^ f ) ~ ( I - r I 73 



4nitv(0)Vft P , V \ I + C- £ Y \ \ \ iu n - £ + £ Y 

E f ( K> . 1 r _ . c afc . 1 ) + ( (4, e _ -L _L_ do 



4vrto(0)^/3 V\l a, iuJ n +£-C Y 17 \\ ^ tuj n + £ - £ Y e k + £ - i v ° 



Q>k<TL>J n 



cL„,dA)-(\ , ' : I „ -4 — ^— 4,4. M )• (bid 



The first and third terms can be combined, and recalling Eq. d52l ) they can be shown to cancel, 



1 1 1 



47Ti ti/(0)Vn j8 V" d - c o fc „ v " ; «v 2 7re^(0) 

since in the steady state we must satisfy I\ = Therefore, this results in: 

= 7- E \ r d 4(W Tr + r TT— A) (I - ^7X7 7X^4. *r 

atlr " M I za; " + £ - £y £k + £ - ir] I / \ I iW„ - £ + £ Y t k - £ + IT) 



(B12) 



(B13) 

I 

In the next step, we will make use of the identity the set of eigenstates of H by {|n)}. With this, one obtains 

de e± L lv ° = ±%n °> (Bi4) r & h +L- °i*> 

f°° 1 

and we briefly outline its proof. It is convenient to use the = / de (m\0\n)de. (B15) 

spectral representation for this purpose and thus enumerate J— 00 

e ± (e m - e„) =F ir l 



14 



After shifting the integration variable e — s- e ± (e m - e„), we 
can further simplify the above expression to yield 



dt- 



1 



P.V. 



[m\0\n) 



(B16) 



de 
e 



±i7T (m|C|n) 



Noting that the principal value integral vanishes concludes the 
proof of Eq. (IB141 >. With this identity in place, Eq. dB 131 ) can 
be brought into its final form 



l_ y l_ 



4. 



C-C Y 
-dl,d r 



-(In 



iuj n — C + Cy a ' a 

R 12 G d„dl (™n) = Q d a dl ( T = °)' 







(B17) 



This demonstrates that the Schwinger-Keldysh approach and 
the effective equilibrium formulation reproduce the same ex- 
pression for the occupancy on the quantum dot. We note that 
the general expressions (53) and (B7) are particularly useful 
since they can be directly applied through numerical proce- 
dures such as numerical RG^^d. 



Appendix C: Derivation of n d using the Keldysh approach 

The basic equation which will constitute the starting point 
of our proof is given by Eq.(15) in Ref.dH The current in lead 
a = ±1 for spin projection a is given by 

x [G< a {t,t 1 ) + fa{e)Gl a {t,t l )} }, 



(CI) 



where we have restricted ourselves to a single level on the dot 
having a spin index a. The occupation number(n £ j (T ) of the dot 
in the state a obeys the following differential equation: 



dn da (t) 
dt 



= — al a0 

— p * — ' 



In steady state (i.e., t=0) we get 

dn da (t) 



dt 



(C2) 



(C3) 



t=o 



This implies 



^/l d<1 /| Im {^ tirQ(Ml) 



[G<(0 l * 1 ) + /.WOS w (0 I ti)]| = 

(C4) 



Since in steady state there exists time translational invari- 
ance, we have G< a (0,h) = G< a (0 - h) and G r aa (0,h) = 
G r a(J (Q — ti). Further we can set r\ — s- and to —> — oo, in the 
sense defined by the open system limit, right from the begin- 
ning. This in turn implies that r Q (0,ti) = T a = 7w(0)i 2 = 
r/2. Let us focus on the first term in the expression above 

A=J2 fjt, J l^e-^TC^Hr)} 

-e^ [GU-ti)]* 
Using the fact [G<,(-h)]' = -G< (ti) 



= --TG< a (0) 

= ir(dt(o)d(o)) = irn d(r . 

Similarly for the second term we get 



(C5) 



B=12 f_Jh J ^lm{e-^f a (e)r a G r aa (-h)} 

a J K J — oo 



= ~T J de/ eff (e,$)^(e). 
Combining the two terms and putting in Eq. (IC4b we get 



(C6) 



n da = I der ff ( e ,$)A^( e ), 



(C7) 



which implies that 

n d = J2 = f def eB (e, $) £ A da (e) 

a •> a 

= J def ff (e^)A d (e). (C8) 

This is identical to the expression obtained via the spectral 
representation in Appendix iBl 
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